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The adhesion of cells is mediated by receptors and ligands anchored in apposing membranes. A cen- 
tral question is how to characterize the binding affinity of these membrane-anchored molecules. For soluble 
molecules, the binding affinity is typically quantified by the binding equilibrium constant Kid in the linear rela- 
tion [RL] — KiD [R] [L] between the volume concentration [RL] of bound complexes and the volume concen- 
trations [R] and [L] of unbound molecules. For membrane-anchored molecules, it is often assumed by analogy 
that the area concentration of bound complexes [RL] is proportional to the product [R] [L] of the area concen- 
trations for the unbound receptor and ligand molecules. We show here (i) that this analogy is only valid for 
two planar membranes immobilized on rigid surfaces, and (ii) that the thermal roughness of flexible membranes 
leads to cooperative binding of receptors and ligands. In the case of flexible membranes, the area concentration 
[RL] of receptor-ligand bonds is proportional to for typical lengths and concentrations of receptors 

and ligands in cell adhesion zones. The cooperative binding helps to understand why different experimental 
methods for measuring the binding affinity of membrane-anchored molecules have led to values differing by 
several orders of magnitude. 



I. INTRODUCTION 

Cell adhesion processes are essential for the distinction of 
self and foreign in immune responses, the formation of tis- 
sues, or the signal transduction across the synaptic cleft of 
neurons 1 1 1. The adhesion processes are mediated by the spe- 
cific binding of receptor and ligand proteins anchored in the 
cell membranes. Because of the importance of these pro- 
cesses, the binding of cells to other cells or to supported 
lipid membranes with anchored ligand molecules has been 
studied intensively with a vaiiety of experimental methods 
m 11 S m. In addition, theoretical models R) E |9l 
[TOl [n [ni [13] [H [Bl [m and experiments on Upid vesi- 
cles with membrane-anchored receptor and ligand molecules 
ifTTl [TSl [191 aim to mimic and capture the specific membrane 
binding processes leading to cell adhesion. 

A central question is how to characterize and measure the 
binding affinity of the membrane-anchored receptor and lig- 
and molecules that are involved in cell adhesion. The binding 
affinity of soluble receptor and ligand molecules can be char- 
acterized by the binding equilibrium constant K^d, defined by 



[RL]iu - i^3D[i?]3D[i]3D 



(1) 



where [i?L]3o is the volume concentration of bound receptor- 
ligand complexes, and [R]3u and [L]3d are the volume con- 
centrations of unbound receptors and unbound ligands in the 
solution. The equilibrium constant K^u is determined by the 
binding free energy of the complex and can be measured with 
standard experimental methods EOl [2ni22]| . An often consid- 
ered two-dimensional analogue for membrane receptors and 
ligands is the quantity 



Kr 



[RL] 
[R][L] 



(2) 



where [RL], [R], and [L] now are the area concentrations of 
bound receptor-ligand complexes, unbound receptors, and un- 
bound ligands |[23ll24ll25l . Since the dimensions of K2D and 



Ki,D are area and volume, respectively. Bell and coworkers 
Il7l [26l suggested that K2D can be estimated as K^d ■ Ic where 
Ic is a suitably chosen 'confinement length', and K^r, is the 
equilibrium constant of soluble counterparts of the membrane 
receptors and ligands obtained by cleaving the membrane an- 
chors. The binding affinity of membrane -anchored receptors 
and ligands have also been directly investigated with fluo- 
rescence recovery after photobleaching Il24l [27l [28l [291 l30l 
and with several 'mechanical methods' |f24l involving mi- 
cropipettes ll25l[3n[32l . hydrodynamic flow chambers |[2l [33l . 
the surface force apparatus |34|, or the biomembrane force 
probe |[35l [36l . However, as pointed out by Dustin and 
coworkers 1241 . the K2D values obtained from fluorescence 
recovery after photobleaching differ by several orders of mag- 
nitude from values measured with mechanical methods. 

Quantifying the affinity of membrane-anchored receptor 
and ligand molecules is complicated by the fact that the bind- 
ing process depends on the local separation and, thus, the 
conformations of the two apposing membranes. We consider 
here a statistical-mechanical model of membrane adhesion in 
which the membranes are described as discretized elastic sur- 
faces and the adhesion proteins as individual molecules dif- 
fusing on these surfaces. We find that thermal shape fluctua- 
tions of the elastic membranes lead to cooperative binding of 
receptors and ligands. The relevant thermal fluctuations occur 
on length scales up to the average separation of the receptor- 
ligand bonds, which is around 100 nm for typical bond con- 
centrations in cell adhesion zones |3|. On these length scales, 
the shape fluctuations are dominated by the bending rigidity 
of the membranes. In our model, the binding cooperativity 
leads to the quadratic dependence 



[RL]^c{K/kBT)llK^,[R]^[L]' 



(3) 



of the bond concentration [RL] on the area concentrations [R] 
and [L] of free receptors and ligands. Here, c ~ 13 is a dimen- 
sionless prefactor, k = kiK2/{ki + K2) is the effective bend- 
ing rigidity of the two apposing membranes with rigidities ki 
and K2, ksT is Boltzmann's constant times temperature, /„e 



2 



the interaction range of the receptor-Hgand bonds, and Kp\ is 
the two-dimensional equiUbrium constant in the case of two 
apposing planar, supported membranes within binding sepa- 
ration of the receptor-ligand bonds. The binding cooperativity 
results from the fact that an increase in the bond concentration 
[RL] 'smoothens out' the membranes, which facilitates the 
binding of additional receptor and ligand molecules. Eq. (|3]l 
holds for typical concentrations and lengths of receptors and 
ligands in cell adhesion zones. The equation implies that K2D 
defined in eq. Q is not constant but depends on the receptor 
and ligand concentration, which helps to understand why flu- 
orescence recovery and mechanical methods to measure this 
quantity can lead to significantly different results, see section 

EH 



11. MODELING BIOMEMBRANE ADHESION 
A. Membrane elasticity 

In our model, the membranes are described as discretized 
elastic surfaces. Each patch of the discretized membranes 
has an area and can contain one receptor or one ligand 
molecule f37l. Molecular models of membranes indicate that 
the whole spectrum of bending deformations is captured if the 
linear patch size a is around 5 nm [381. 

The elasticity of lipid membranes in general depends on the 
bending rigidity k ||39l and the tension a. An important length 
scale is the 'crossover length' y/ n/a ll40l . The membrane ten- 
sion dominates over the bending energy on length scales larger 
than nja, while the bending energy dominates on smaller 
length scales. For typical values of the bending rigidity k of 
lipid bilayers around 10~^^ J ~ 25 fcsT [41 1 and tensions of 
a few [iSIm ll42ll . the crossover length \fnja attains values 
of several hundred nanometers. In addition, the elasticity of 
cell membranes is affected by the actin cytoskeleton on length 
scales larger than the distance of the cytoskeletal membrane 
anchors, which is around 100 nm 1 1 1. 

In cell or membrane adhesion zones, the relevant shape 
fluctuations occur on length scales up to the average distance 
of the receptor-ligand bonds, since the bonds locally constrain 
the two membranes. Typical values for the average bond dis- 
tance in cefl adhesion zones range from 50 to 100 nm 13], 
which is significantly smaller than the crossover length n/a 
estimated above, and smaller than or equal to the distance of 
the cytoskeletal membrane anchors. The relevant shape fluctu- 
ations in the adhesion zones on length scales up to the average 
bond distance are therefore dominated by the bending energy. 
In our model, the bending energy has the form l(37l 



(4) 



where li is the local separation of the apposing membrane 
patches i in the adhesion zone. Ad is the discretized Lapla- 
cian given in appendix A, and k = kiK2/(ki + K2) is the 
effective bending rigidity of the two membranes with rigidi- 
ties Ki and K2- If one of the membranes, e.g. membrane 2, 



is a planar supported membrane, the effective bending rigid- 
ity n equals the rigidity ki of the apposing membrane since 
K = niK2/{ni + K2) ~ 1^1 for large K2. 



B. Receptor-ligand interactions 

We consider here a membrane with a single type of receptor 
molecules apposing a second membrane with complementary 
ligands, see fig.[TJa). A receptor can bind to a ligand molecule 
in our model (i) if the ligand is located in the membrane patch 
apposing the receptor, and (ii) if the membrane separation U is 
close to the length lo of the receptor-ligand bond. For simplic- 
ity, we describe the receptor-ligand interaction by the square- 
well potential 

V{h) = -U for lo - Ce/2 <k<lo + Zwe/2 

= otherwise (5) 

where C/ > is the binding energy and l^^ the binding range 
of a receptor-ligand complex. The binding range is the dif- 
ference between the smallest and the largest local membrane 
separation at which the molecules can bind. The binding 
range l^^ depends on the interaction range of the two bind- 
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FIG. 1: (a) A membrane with receptor molecules (top) interacting 
with ligands embedded in an apposing membrane (bottom). A re- 
ceptor can bind a ligand molecule if the local separation of the mem- 
branes is close to the length of the receptor-ligand complex. - (b) In 
our model, the attractive interactions between the receptor and lig- 
and molecules lead to an effective single-well adhesion potential Kf 
of the membranes. The depth U^f of the effective membrane binding 
well depends on the concentrations and binding affinity of receptors 
and ligands, see eq. l|7j. The width i„c of the binding well is equal to 
the binding range of the receptor-ligand interaction, see eq. |5|. 
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ing sites on the receptor and ligand molecules, on the flexibil- 
ity of these molecules, and on the flexibility of the membrane 
anchoring. For the rather rigid protein receptors and ligands 
that typically mediate cell adhesion, the binding range may be 
around 1 nanometer 

The interaction energy of the membranes is then described 

by 

Hint{l, n, m} = ^ n^miVik) (6) 

i 

Here, rii = 1 or indicates whether a receptor is present 
or absent in membrane patch i, and = 1 or indicates 
whether a ligand is present or absent in patch i of the appos- 
ing membrane. The configurational energy of the membranes 
H{1, n, to} = 'Hc\{l} + 'Hint{l, n, to} is the sum of the elastic 
energy Q and the interaction energy (j6]l. 

III. EFFECTIVE ADHESION POTENTIAL OF THE 
MEMBRANES 

The equilibrium properties of our model can be derived 
from the free energy !F = — fc^TlnZ, where Z is the par- 
tition function, /c^ is Boltzmann's constant, and T is the tem- 
perature. The partition function Z is the sum over all possible 
membrane configurations, with each configuration {l,n,m} 
weighted by the Boltzmann factor exp [— 7Y{Z, n, to}/A:bT]. 
In our model, the partial summation in Z over all possible 
distributions to and n of receptors and ligands can be per- 
formed exactly, which leads to an effective adhesion poten- 
tial, see appendix A. The effective adhesion potential of the 
membranes is again a square-well potential of the form Q, 
with the same binding range l^^f. as the receptor-ligand inter- 
action, but with an effective potential depth Uef that depends 
on the concentrations and binding energy U of receptors and 
ligands, see fig.[TJb). The concentrations of receptors and lig- 
ands in biological or biomimetic membranes are several or- 
ders of magnitude smaller than the maximum concentration 
1/a^ ~ 4 • W^/iim^ in our discretized membranes with patch 
size a ~ 5 nm. For these small concentrations, the effective 
potential depth is 

C/ef «fcBr[i?][L]a2e^/'=«^ (7) 

as shown in the appendix. The equilibrium behavior of our 
model thus can be determined from considering two lipid 
membranes interacting via an effective adhesion potential 
with well depth Ucf. 

IV. MEMBRANE FRACTION WITHIN BINDING RANGE 
OF RECEPTORS AND LIGANDS 

A receptor molecule can only bind an apposing ligand 
molecule if the local membrane separation is comparable to 
the length of the receptor-ligand complex. A central quantity 
in our model therefore is the fraction of apposing mem- 
brane patches i with a separation li within the binding range 



lo ± ?we/2 of the receptor-ligand interaction ([s]). Our goal is 
to determine the equilibrium concentration [RL] of receptor- 
ligand bonds. We find that this concentration is proportional 
to Pb, and proportional to the concentrations [R] and [L] of 
unbound receptor and ligand molecules, and given by 

[RL]^Pb[R][L]a^e^^''^^ , (8) 

see appendix. The quantity a^e'^^'^^^ in this equation can be 
understood from considering first two planar, parallel mem- 
branes supported on rigid substrates, e.g. two supported mem- 
branes in the surface force apparatus ||34ll43l . If the separation 
of the two membranes is close to the length lo of the receptor- 
ligand complex, we have Pb = 1. A comparison with eq. ([T]| 




FIG. 2: The roughness of the membranes decreases with increas- 
ing concentration of receptor-ligand bonds (top to bottom) since 
the bonds constrain membrane shape fluctuations. The fraction Pb 
of membrane segments within binding separation of the receptors 
and ligands therefore increases with the bond concentration (see 
eq. l|14[l), which leads to binding cooperativity: The binding of re- 
ceptors and ligands 'smoothens out' the membranes and, thus, facil- 
itates the formation of additional receptor-ligand bonds. For clarity, 
unbound receptor and ligand molecules are omitted in the cartoons. 
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then indicates that the quantity 



(9) 



can be interpreted as the two-dimensional binding equilibrium 
constant of the receptors and ligands in the case of planar 
membranes with separation lo- If we now use the definition 
(|2]i for K2D, we obtain 



K2D = PbK^ 



(10) 



However, thermal shape fluctuations of flexible membranes 
not supported on rigid substrates can lead to values of Pi, 
much smaller than 1 . We will show here that Pi, depends on 
the concentrations of the receptors and ligands, which results 
in cooperative binding, see fig. |2] As an equilibrium quantity, 
the area fraction Pi, of the membranes within receptor-ligand 
binding range is determined by the effective adhesion poten- 
tial shown in fig. [ijb). 



V. LONG RECEPTOR-LIGAND COMPLEXES 

We consider first the case in which the length of the 
receptor-ligand complexes is larger than the thermal mem- 
brane roughness. The fluctuating membranes then do not 
touch each other, and the repulsive hard-wall interaction of 
the lipid membranes can be neglected. We will show in the 
next section that this case applies to biological receptor-ligand 
complexes with typical lengths between 15 and 40 nm ll44l . 
which are much larger than the binding range of the re- 
ceptor and ligand molecules. Scaling arguments indicate that 
the central parameter that affects Pi, is the rescaled effective 
potential depth |45 1 



U = UefKll^/{kBTy 



(11) 



in the case of long receptor-ligand complexes. The fraction 
Pb — Pb{u) of membrane patches with a separation within 
binding range of the receptors and ligands thus is a function 
of a single parameter u in this case. 



A. Linear regime 

For relatively small values of the rescaled effective potential 
depth u, the membrane fraction Pi, within binding range of the 
receptors and ligands is linear in u and behaves as 



Pb cu = ckII^ Uef/ikBTy 



(12) 



The numerical prefactor c — 13 ± 1 in this equation can be 
determined from Monte Carlo simulations, see figs. |3] and |4] 
With eq. (|7]), we obtain 



Pb~ciK/kBT)llMR][L] 



(13) 



and, thus, a linear dependence of Pb on the concentrations [R] 
and [L] of unbound receptors and ligands. Our central result 
Q then follows directly from inserting eq. ( 13 1 into eq. (|8]l. 
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FIG. 3: Area fraction Ph of the membranes within binding range of 
receptors and ligands in the case of long receptor-ligand complexes 
in which the hard-wall repulsion of the lipid membranes is negligi- 
ble, see text. Here, u = Ua ^,11,^/ {ksTY is the rescaled effective 
potential depth, which is the central parameter governing Pj, in this 
case. The Monte Carlo data points are well fitted by the function 
Pi, ~ (?i + C2V? + Csu^)/ {ci + u + C2U^ + CiU^) with the three 
fit parameters ci ~ 0.070, C2 — —0.32, and C3 ~ 0.50, see full 
line. - (Inset) For < Pb < 0.6, the Monte Carlo data points 
can be described by Pt — u/{ci + u) with the single fit parame- 
ter ci ~ 0.071. This single-parameter function is also a reasonable 
approximation over the whole range of Pi, values, see dashed line. 
The linear relation l |12[ > results from an expansion of this function for 
small values of u. 



The linear behavior ( [T2| ) does, in fact, hold for a wide, bio- 
logically relevant range of concentrations and bending rigidi- 
ties. To see this, we consider the relation 



Pb « Vc{i^/kBT)lURL] 



(14) 



obtained from eq. ( [T2| ) and Uef — kBT[RL]/Pb, which fol- 
lows from eqs. (j7]i and ([8]l. The linear relation ( [T2) i is valid 
for small Pb < 0.2, see fig. [3] Typical bond concentrations in 
cell adhesion zones are around [RL] — 100//im-^ |3|, while 
the binding range of receptor and ligand proteins can be esti- 
mated as Z„e = 1 nm. Typical values for the bending rigidities 
of lipid bilayers are around 25 /csT PTl , which implies an 
effective rigidity k of 12.5kBT for two apposing membranes, 
see text below eq. (|4]). Because of its embedded and attached 
proteins, the bending rigidities of biological membranes may 
be larger, e.g. by a factor 2. For the effective rigidity of lipid 
bilayers, we obtain the estimate Pb ~ 0.13 from eq. ^T4[ , and 
for a 2-fold increased effective rigidity, we obtain Pb ~ 0.19. 
Both estimates are within the range of Pb values for which the 
linear relation ( [T2j i is valid. 



B. Nonlinear regime 

Our model is not limited to the linear regime considered in 
the previous section. The full functional dependence of Pb on 
the rescaled effective potential depth u defined in eq. 
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be determined from Monte Carlo simulations, see fig. |3] The 
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membrane fraction Pi, within receptor-ligand binding range is 
linear in u for small values of u, and increases to 1 for large 
values of u. We find that the Monte Carlo data can be fitted 
well by 



Ph 



Cl 



(15) 



in the range < Pf, < 0.6 with the single fit parameter ci = 
0.071 ± 0.002, see inset of fig. |3] In addition, eq. ([T5]l is a 
reasonable approximation over the whole range of Pi, values, 
see dashed line in fig. [3] Inserting eq. ( fTS) into eq. ([8]l leads to 



[RL] 



(16) 



This equation generalizes eq. Q beyond the linear approxi- 
mation valid for Pb < 0.2. 



VI. SHORT RECEPTOR-LIGAND COMPLEXES 

Lipid vesicles with anchored receptor molecules are im- 
portant biomimetic systems for cell adhesion. In these sys- 
tems, the receptors either bind to ligands anchored in other 
vesicles ifTSl . or to ligands anchored in supported membranes 
ifTTl [191. In principle, the length of the receptor and ligand 
molecules in these systems can be varied. For short receptor- 
ligand complexes, the hard-wall interaction of the lipid mem- 
branes is an important aspect of adhesion. The hard-wall in- 
teraction leads to an entropic, fluctuation-induced repulsion of 
the membranes |46|, and to an unbinding transition at a finite 
strength of the attractive interactions between the receptors 
and ligands ll47l . In our model, the strength of these attractive 
interactions is captured by the depth Uef of the effective ad- 
hesion potential shown in fig.[l] In the interplay between en- 
tropic repulsion and attractive receptor-ligand interactions, the 
membranes will be bound for potential depths Ugf > Uc where 
Uc is the critical interaction strength of the binding transition. 
The membranes are unbound for potential depths J/gf < Uc- 

The critical potential depth Uc depends (i) on the width Z„e 
of the potential well, which is identical with the binding range 
of the receptor-ligand interaction (|5]l, and (ii) on the separation 
li of the potential well from the hard wall, see fig.[TJb). For 
h > 'we, we have previously obtained the relation 



Uc 



b{kBTY 



(17) 



1 ^we 



from scaling arguments and Monte Carlo simulations, with the 
numerical prefactor b — 0.025 ± 0.002 ll45l . With increasing 
li, the critical potential depth Uc goes to zero since the effect 
of the entropic repulsion decreases. 

In the bound state of the membranes with U^f > Uc, the 
membrane fraction P^ within binding range of the receptors 
and ligands increases with increasing values of U^f, see fig. [4] 
As in the previous section, we consider here first the linear 
regime for values of U^i close to Uc- We find that eq. ( [T2| , 
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FIG. 4: Area fraction Pt of the membranes within binding range of 
receptors and ligands as a function of the rescaled effective potential 
depth It for different values of lv,c/li- Here, Zwc is the width of the 
effective potential well shown in fig.[TJb), and h is the separation of 
this well from the 'hard wall' atl — where the two lipid membranes 
are in contact. The Monte Carlo data for iwe/'i ~ are the data 
shown in fig. [3] The full line interpolating these data is the three- 
parameter fit function given in the caption of fig. [3] The dashed lines 
are quadratic fits to extrapolate the data for iwc/'i ~ 0.5, 1 and 2 to 
Pi, = 0. From these fits, we obtain the following estimates for the 
slope Pfc = dPb/du al Pt, = 0:22 ±2 for L^/h = 0.5, 24 ± 2 for 
Ivic/h = 1, and 19 ± 2 for l^s/h = 2. From a similar quadratic fit of 
the data for Isxe/h ~ 0, we obtain the slope P5(0) = 13 ± 1, which 
is our estimate for the parameter c in eq. \12\ - From the three values 
of Pj for — 0, 0.5 and 1 at Pb = 0, we obtain the estimate 

d = 20 ± 5 for the parameter d in eqs. |T9| and l |20[ (. 



which holds for large values of h, can be generalized to 

Pb = for f/ef < Uc 

« Jl + g(Ce/'l)) (f/ef - Uc)/{kBTf for C/ef > Uc 

(18) 

with the same numerical prefactor c = 13 ± 1 as eq. ( [T2| , and 
a function .g(Ce/'i) that goes to zero for small Ce/'i- In the 
limit of large li, eq. ( [T8| is then identical with eq. ( [T2) l. In this 
limit, the function g{lwe/h) can be approximated by 



(19) 



which can be unterstood from a Taylor expansion of g{x) 
around x = 0- From the Monte Carlo data shown in fig. |4] 
we obtain the value d = 20 ± 5 for the numerical pref- 
actor in eq. ( [T9) l. For [/ef > Uc, which is equivalent to 
ii'pi[i?][P] > Uc/ksT, eq. ^ thus generaUzes to 

[RL] « c{K/kBT)ll^{l + dl^,/h)x 

X AV[P][P](i^pi[i?][P] " Uc/keT) (20) 



For /<'pi[i?][L] < Uc/ksT with Uc given in eq. ( 17 1, we have 
[RL] = since the membrane fraction Pi, within binding 
range of the receptors and ligands vanishes in this case. 

The lengths lo of receptor-ligand complexes in cell con- 
tact zones range from 15 to 40 nanometers |44|. For these 
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receptor-ligand complexes, the separation h ~ lo — Ivie/'^ of 
the effective potential well in fig.[TJb) from the 'hard wall' at 
I = is much larger than the width /^e of the well, which 
can be estimated to be of the order of 1 nm, see section [nB] 
In this case, relation (20 1 is practically identical with relation 
(|3]l, i.e. the two relations are identical within the numerical er- 
rors of the parameters b, c and d at these values of li and /^ve> 
which implies that the hard-wall interaction is negligible for 
biological receptor-ligand complexes with a length between 
15 and 40 nm. 

This conclusion can be confirmed by considering scaling 
estimates for the membrane roughness = (If) — (k)^. 
Here, (k) and (If) are statistical averages of the local mem- 
brane separation 1^ and its square If. For l^g <C h, the rough- 
ness is approximately ~ 0A4:^/k^^TJ(K\R^) BSl 
and, thus, attains the value 4 nm for the typical bond con- 
centration [RL] ~ 100//im^ |3| and the effective rigidity 
K — kiK2/{ki + K2) = 12.5kBT of two apposing mem- 
branes with the bending rigidities ki — K2 = 25A;bT. This 
estimate for the roughness is much smaller than the length 
of the receptor-ligand complexes. The steric repulsion of the 
membranes arising from direct membrane-membrane contacts 
therefore is negligible, and the average separation (k) of the 
membranes is equal to the length lo of the receptor-ligand 
complexes. 



VII. DISCUSSION AND CONCLUSIONS 

We have shown here that a central quantity in cell adhe- 
sion is the fraction of the apposing membranes within bind- 
ing range of receptors and ligands. In equilibrium, the aver- 
age separation of two membrane segments bound by receptor- 
ligand complexes is equal to the length of the complexes, pro- 
vided that the steric repulsion of the lipid membranes and 
other repulsive interactions, e.g. from large glycoproteins, are 
negligible. The fraction of membrane patches with a sep- 
aration within the receptor-ligand binding range then depends 
on the thermal roughness of the membranes, which in turn is 
affected by the concentrations of the receptors and ligands. 
For small concentrations of receptors and ligands, the fraction 
Pfc of the membranes within binding range is proportional to 
[R][L], see eq. {13}. Since we have [RL] = Pb[R][L]Kpi ac- 
cording to eq. (jsjl, this proportionality leads to the quadratic 
dependence ^ of the bond concentration [RL] on the concen- 
trations [R] and [L] of free receptors and ligands, which indi- 
cates cooperative binding. The linear relation ( [T3| ) between 
Pb and is valid for small Pb < 0.2. The more general 

relation ( [TSj ), which is a good approximation over the whole 
range of Pb values, leads to eq. ( [T6] l. 

Our results may help to understand why experimental val- 
ues for K2D defined in eq. Q obtained with the fluorescence 
recovery method are several orders of magnitude larger than 
the values obtained with the micropipette method 1241 . In flu- 
orescence recovery experiments, K2D is measured in the equi- 
librated contact zone of a cell adhering to a supported mem- 
brane with fluorescently labeled ligands. In micropipette ex- 
periments, in contrast, K2D is measured for initial contacts be- 



tween two cells. Dustin and coworkers 112411 have pointed out 
that the different orders of magnitude of K2D obtained with 
these two methods can be partly understood from different 
contact areas. In the micropipette experiments, large mem- 
brane protrusions such as microvilli may lead to an actual con- 
tact area Ac that is only a few percent of the observed contact 
area. However, even correcting for a significantly smaller ac- 
tual contact area Ac in the micropipette experiments still leads 
to values of K2D that are 3 to 4 orders of magnitude smaller 
than the K2D values from fluorescence recovery experiments, 
see fig. 2 in ref. |24J. 



We suggest that this orders-of-magnitude gap can be fur- 
ther closed by considering PbAc as the relevant quantity for 
receptor-ligand binding, rather than Ac- The quantity PbAc 
is the fraction of the actual contact area Ac in which the 
two membranes are within binding separation of the receptor- 
ligand bonds. Since the number of receptor-ligand bonds is 
proportional to PbAc, differences in this quantity in experi- 
mental setups translate directly into differences in K2D, see 



eq. (10 1. There are two significant differences between the flu- 
orescence recovery and the micropipette experiments. First, in 
the equilibrated contact zones of the fluorescence recovery ex- 
periments, the bond concentration [RL] is enriched, compared 
to the bond concentration for initial cell contacts probed in the 
micropipette experiments. This enrichment results from a dif- 
fusion of free receptor and ligand molecules into the contact 
zone, in which the molecules can bind. According to eq. ( [T4| , 
an increase in [RL] by a factor of 100, which is not unrealis- 
tic Q, leads to an increase in Pb by a factor of 10, and thus 
explains one order of magnitude in the observed difference of 
K2D values from fluorescence recovery and micropipette ex- 
periments. Second, the average separation of cell membrane 
and supported membrane in the fluorescence recovery experi- 
ments is close to the length of the receptor-ligand complexes. 
For the initial cell-cell contacts in the micropipette experi- 
ments, in contrast, the average separation of the membranes 
will deviate from the length of the receptor-ligand complexes, 
e.g. because of large glycoproteins that eventually will diffuse 
out of the contact zone. The effect of this deviation of the 
average membrane separation on Pb is more difficult to as- 
sess, but may easily account for an additional 1 or 2 orders 
of magnitude difference in PbAc between fluorescence recov- 
ery and micropipette experiments. The membrane fraction Pb 
within binding range depends sensitively (i) on the thermal 
membrane roughness, and (ii) on the difference between the 
average membrane separation and the length of the receptor- 
ligand complexes. Pb is small if the difference between av- 
erage membrane separation and complex length exceeds the 
thermal membrane roughness, which also leads to small val- 
ues of K2D, see eq. ( [T0| . 

We have made several simplifying assumptions in our 
model. One of these simplifications is the square-well form 
(|5]l for the receptor-ligand interaction. A convenient aspect of 
the square-well interaction is that the effective adhesion po- 
tential of the membranes, which results from an integration 
over the receptor and ligand degrees of freedom in the parti- 
tion function, has the same square-well form, with an effective 
depth Usf that depends on the concentrations and on the equi- 



7 



librium constant Kpi of the receptors and ligands in the case 
of planar membranes, see eqs. (|7]i and However, the effec- 
tive potential can also be calculated for other functional forms 
of the receptor-ligand interaction, e.g. for a Gaussian form. In 
general, two important parameters of the receptor-ligand inter- 
action are the width Zwe and depth U of the potential well. The 
width l^f. is affected by the interaction range of the binding 
sites on the receptors and ligands, by the flexibility of these 
molecules, and by the membrane anchoring, while the depth 
U is directly related to the equilibrium constant Kp\ for pla- 
nar membranes, see eq. (j9]l. We expect that the square-well 
potential (j5]l and other functional forms of the receptor-ligand 
interaction lead to rather similar results for comparable values 
of the width and depth of the potential wells. 



We have argued in section II A that the elasticity of cell 
membranes is dominated by their bending rigidity on the 
length scales up to 100 nm relevant here. For typical tensions 
a of the membranes, the crossover length k/ct, above which 
the tension dominates over the bending energy, is clearly 



larger, see section II A However, these length scales are only 



slightly smaller or comparable to the average separation of 
the cytoskeletal anchors in cell membranes. On the one hand, 
the anchoring to the cytoskeleton may suppress thermal fluc- 
tuations on length scales larger than the average separation 
of the anchors. On the other hand, active processes within 
the cytoskeleton may increase membrane shape fluctuations 
||49l[50ll . In general, active cell processes and inhomogeneities 
may perturb the homogeneous equilibrium situation consid- 
ered here. However, our results are still applicable to mem- 
brane regions in which the concentrations of bound and un- 
bound receptors and ligands are locally equilibrated. In prin- 
ciple, vesicles and supported membranes with anchored re- 
ceptors and ligands are excellent model systems to test our 
theoretical results without the complications of the cell cy- 
toskeleton and of active biological processes. 



APPENDIX A: EFFECTIVE ADHESION POTENTIAL 

In this appendix, we derive the effective potential shown in 
fig. [ijb). We perform the calculations in the grand-canonical 
ensemble in which the concentrations of receptors and ligands 
are adjusted by the chemical potentials and /i^. The chem- 
ical potentials are free energy differences between a patch of 
size of our discretized membranes that contains a receptor 
or ligand molecule and a membrane patch without receptor or 
ligand. 

In the grand-canonical ensemble, the configurational en- 
ergy of the membranes is 371 



H{1, n, m} = Hci{l} + Hint{^ n, m} - ^ (n^^^j + rriifiL) 

(AI) 

with the elastic energy Q and interaction energy (j6|l. The 
elastic energy (j4|l depends on the mean curvature {Adli)/a'^ of 
the separation field k, with the discretized Laplacian A^^^ = 
ki + U2 + hi + hi — 4:1 i- Here In to li4 are the membrane 
separations at the four nearest-neighbor patches of membrane 
patch i on the quadratic array of patches. 



The equilibrium properties of our model can be determined 
from the free energy JF = —ksThiZ, where Z is the parti- 
tion function, fc^ is Boltzmann's constant, and T is the tem- 
perature. The partition function 



n 



HE 



HE 



-H{l,n,m}/kBT 



(A2) 



is the sum over all possible membrane configurations, with 
each configuration {I, n, m} weighted by the Boltzmann fac- 
tor exp [— 7i{/, n, mj/fesT]. In our model, the summations 
in Z over all possible distributions m and n of receptors and 
ligands can be performed exactly, which leads to ifSl lJTl 



n 



dk 



,-(Hci{(}+a^ T.i V,i(li))/kBT 



(A3) 



with the effective adhesion potential 

V,,{k) = - ^In [Co + (e-^('»)/'=-^ - l)e(A'H+ML)/fc«T 

(A4) 

and Co = (1 + e^«/'=B^) (l + e'^^/''^'^). For the receptor- 
ligand interaction (|5]l, the effective adhesion potential is a 
square-well potential with the same width /wc and the depth 



C/ef = 



■In 



(A5) 

The total concentrations of receptors and ligands follow 
from partial derivates of the free energy with respect to the 
chemical potentials: 



[R] + [RL] 
[L] + [RL] . 



1 



1 
1 



(A6) 
(A7) 



Here, A denotes the membrane area. The concentration [RL] 
of receptor-ligand bonds is obtained from a partial derivative 
with respect to the binding energy U of the bonds: 



1 dF 
AdU 



[RL]: 

These three equations lead to 



1 



(A8) 



[R] = 



1 



(1 - Pb) e^«/'=^'^ 
1 + e'^^/'^^'^ 



\ -)_ gPij/fcsT _|_ gfj.L/kBT _j_ ^{U+fiR+fj.L)/kBT 



(A9) 



[L] = 



(1 - Pb) e 



f^h/kBT 



1 + ef^i'/^'BT 



1 _|_ (.P-R/ksT _|_ gllL/kBT _|_ f^(U+flR+tlL)/kBT 



(AlO) 
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and 



(All) 

where Pi, is equilibrium fraction of membrane patches with 
a separation li within the binding range lo — lv,e/'^ < U < 
lo + ^we/2 of the receptor-ligand interaction (jsji. 

The typical concentrations of receptors and ligands in cell 
membranes up to several hundred molecules per square mi- 
cron are significantly smaller than the maximum concentra- 
tion ~ 4.1 • 10''//im-^ in our model. This implies 



in eqs. ( A9 1 to ( Al 1 1. With these relations, we obtain 



[R] « 



P-R/ksT 



[L] « 



iJ-L/kBT 



(A12) 



and 



[RL] 



1 



p g{U + fJ.R+l-lL)/kE 



Pb[R][L] a'e 



2 U/ksT 



(A13) 



The effective potential depth ( |A5| l then simplifies to 

f/ef 



(A14) 



APPENDIX B: MONTE CARLO SIMULATIONS 

The area fraction Pt, of the membrane within the well of 
the effective potential can be determined with Monte Carlo 



simulations |]8] |37l |45]- It is convenient to use the rescaled 
separation field Zi = [li/ a)yj n/ [kBT) in the simulations. 
The configurational energy then has the form ^{{z^/kBT = 
X]Jf (Ad^i)^ + a'^Vfif{zi)/kBT] where Vef is the effective 
potential shown in fig. [ijb). In the Monte Carlo simulations, 
local moves are attempted in which the rescaled separation 
Zi of the membrane patch i is shifted to a new value z,; + ( 
where is a random number between —1 and 1. Follow- 
ing the standard Metropolis criterion [51], a local move is al- 
ways accepted if the change AH in conformational energy is 
negative, and accepted with the probability exp(— ATi/fesT) 
for AH > 0. We perform simulations with up to 5 • 10^ 
attempted local moves per site i and membrane sizes up to 
= 160 X 160 patches. The membrane size is always cho- 
sen to be much larger than the lateral correlation length of the 
membranes. Thermodynamic averages of the fraction Pi, of 
membrane patches bound in the potential well then do not de- 
pend on the finite system size. The Monte Carlo data shown 
in figs. |3] and |4] are from simulations with the rescaled width 
z„e = 1 of the potential well. Further details of our Monte 
Carlo simulations are described in ref. 1371. 
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